On the Mass Profile of Galaxy Cluster CL 0024+1654 Inferred 

from Strong Lensing 



Paul R. Shapiro 
Department of Astronomy, University of Texas, Austin, 78712 

shapiro@astro . as . utexas . edu 

and 

Ilian T. Iliev 

Department of Physics, University of Texas, Austin, 78712 
iliev@astro . as . utexas . edu 

Received ; accepted 



submitted to ApJ Letters 



- 2 - 



ABSTRACT 

Observations of a flat density profile in the cores of dark-matter-dominated 
halos on the two extremes of mass for virialized objects in the universe, dwarf 
galaxies and galaxy clusters, present a serious challenge to the current standard 
theory of structure formation involving Cold Dark Matter (CDM). By contrast, 
N-body simulations of halo formation in the latter indicate density profiles which 
are singular and steeply rising towards the center. A flat-density core on the 
cluster scale is indicated by gravitational lensing observations, most significantly 
by the strong-lensing measurements of CL 0024+1654 by the Hubble Space Tele- 
scope. A recent re-analysis of this cluster has suggested that a uniform-density 
core is not demanded by the data, thereby eliminating a significant piece of 
the conflict between the observations and the CDM theoretical predictions. We 
show here, however, that the singular mass profile which that analysis reports 
as consistent with the lensing measurements of CL 0024+1654 implies a velocity 
dispersion which is much higher than the measured value for this cluster. 

Subject headings: cosmology: gravitational lensing - cosmology: theory - dark 
matter - galaxies: clusters: individual (CL 0024+1654) - galaxies: halos 

1. Introduction 

There has been a lot of recent controversy concerning the density profiles of the dark 
matter halos of virialized cosmological structures, from dwarf galaxy to galaxy cluster scale. 
N-body simulations of the formation of "virialized" dark matter halos associated with 
galaxies and clusters in a Cold Dark Matter model were found to be well-fit by a simple, 
universal form for the variation of mass density p with radial distance r from the center of 
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mass, given by 

PNFw(?") = "T - / TTT"^ / ^> (!) 

(r/r 5 )(l + r/r s y 

where r s is some characteristic radius which separates the two asymptotic power-law slopes, 
p oc r^ 1 at r << rs and p oc r~ 3 at r >> rs, and ps is a characteristic density which 
reflects the mean density of the universe at the epoch of halo formation (Navarro, Frenk, 
& White 1997; "NFW"). More recent N-body simulations of higher resolution obtain halo 
profiles which agree with the NFW profile at large radii but have an even steeper inner 
slope, p oc r~ 15 (Moore et al. 1999). 

This prediction by the standard CDM model of singular density profiles for cosmological 
halos is apparently in conflict with the observed mass distributions inside dark-matter- 
dominated halos on two extremes of the halo mass function, dwarf galaxies and galaxy 
clusters. As a result, these observations and their interpretation have recently come under 
intense scrutiny. In work going back to Flores & Primack (1994), Moore (1994), Burkert 
(1995) and more recently in Kravtsov et al. (1998), Moore et al. (1999), and McGaugh & 
de Blok (1999), these predictions of singular halos were found to be in conflict with the halo 
density profiles derived from the observed rotation curves of dark-matter-dominated dwarf 
galaxies. The latter were found, instead, to be better fit by density profiles with flat cores. 
Recently, the universality of this observational requirement that dark-matter-dominated 
galactic halos possess a uniform- density core has been challenged on the grounds that the 
rotation curves are not generally resolved well enough in the centers to distinguish the 
slowly rising rotation curve which results from a mass profile with a flat core from the more 
rapid rise which results from a cuspy profile (van den Bosch et al. 2000). Initially, this 
challenge appeared to apply primarily to LSB galaxies, leaving the case for uniform- density 
cores in dwarf galaxies still strong. For example, the well-resolved rotation curves of 
the nearby dwarf galaxies DDO 154 and NGC 3109 still demanded a flat density core 
(van den Bosch et al. 2000). Very recently, however, this challenge has been extended to 
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the interpretation of dwarf galaxy rotation curves, as well (van den Bosch & Swaters 2000). 
While the latter are still generally better fit by halo profiles with a flat core, it can no longer 
be stated with much confidence that the observations demand such profiles. Until better 
resolution data become available, profiles with p oc r~ a near the center, with < a < 1, all 
seem to yield reasonable fits to the data. Hence, while it appears that a significant conflict 
remains between the CDM halos predicted by the highest resolution N-body simulations 
to date (which imply p oc r~ L5 at the center) and the dwarf galaxy rotation curves, the 
current observations are unable to discriminate effectively between halos with inner profiles 
as cuspy as the NFW profile and those with a uniform-density core, instead. 

On the galaxy cluster scale, the case for a uniform-density core has been made 
most convincingly using observations of strong gravitational lensing, where the images of 
background galaxies are distorted by the cluster mass to form arcs and multiple images. 
To date, the most spectacular example of such observation is a Hubble Space Telescope 
image of multiple arcs produced by the cluster CL 0024+1654 at z — 0.39. The relaxed 
structure of CL 0024+1654, the absence of a single, central, dominant cluster galaxy, and 
the presence of an easily-identified multiply-imaged background galaxy make this cluster 
a particularly good candidate for a determination of the halo mass profile from lensing 
analysis. According to Tyson, Kochanski, & Dell'Antonio (1998), the observations of this 
cluster require a halo mass profile with a uniform-density core, in strong conflict with the 
predicted cuspy profile of NFW. 

These conflicts between the CDM N-body results and the observations of dwarf 
galaxy and cluster halo profiles have stimulated a vigorous re-examination of the 
theoretical underpinnings of the CDM model, including a number of suggestions for 
a variation of the microscopic properties of CDM which might serve to produce halos 
with uniform- density cores while retaining the more successful aspects of the original 
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CDM model. These include suggestions that the dark matter is nongravitationally 
self- interacting (Spergel & Steinhardt 2000), warm (e.g. Sommer-Larsen & Dolgov 2000; 
Colin, Avila-Reese, & Valenzuela 2000; Hannestad & Scherrer 2000), fluid (Peebles 2000), 
decaying (Cen 2000), repulsive (Goodman 2000), fuzzy (Hu, Barkana, & Gruzinov 2000), 
and annihilating (Kaplinghat, Knox & Turner 2000). In view of the importance of the 
dwarf galaxy and cluster halo profile observations in constraining the theoretical models, it 
is somewhat disappointing that the conclusions based upon dwarf galaxy rotation curves 
are currently so ambiguous with regard to the question of the uniform-density core. This 
makes the conclusion of Tyson et al. (1998) regarding the core in CL 0024+1654 all the 
more critical. 

Recently, a new study of this cluster by Broadhurst et al. (2000) has reached a 
conclusion opposite to that of Tyson et al. (1998) regarding the consistency of the observed 
mass profile with the density profiles predicted by the N-body simulations of cluster 
formation in the standard CDM model. They find that the NFW mass profile is consistent 
with the lensing data. If this is correct, then the case against the standard CDM model is 
significantly weakened. The purpose of this letter is to point out that the fit by Broadhurst 
et al. (2000) of the cluster lensing data with a mass distribution which follows the NFW 
profile implies a cluster velocity dispersion which is much larger then the value measured 
for this cluster by Dressier et al. (1999) of oy = 1150 km s -1 . 



2. Observational Results for the Cluster Mass Profile 

The projected dark matter density profile which Tyson et al. (1998) found by modeling 
the lensing data for CL 0024+1654 within the arcs at radius r arcs ~ 100 h^kpc is well fit by 

K(l + r ] y 2 ) 

Z(V) = (1+ 3,2)2-, . ( 2 ) 
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where y = r/r corc , K = 7900 ± 100 h M pc~ 2 , r corc = 35 ± 3 h' 1 kpc, rj = 0.57 ± 0.02, and h 
is the Hubble constant in units of lOOkms^Mpc -1 . Additionally, Tyson et al. (1998) rule 
out at a great confidence level the possibility of a good fit of the observed mass distribution 
by the NFW profile. 

A recent paper by Broadhurst et al. (2000) obtained, instead, a good fit to the 
lensing data with a total mass distribution given by the NFW profile in equation (1), with 
r s ~ 400 h~ l kpc and 5 C « 8000, where 5 C = ps/p C rit(z), and p cr i t {z) = 3H 2 /(8nG), the 
critical density of the universe at the cluster redshift. This 5 C is directly related to the NFW 
concentration parameter c = r 2 oo/rs, where r 2 oo is the radius within which the average 
density is 200 times this critical density, according to 

_ 200 

° c 3 [m(l + C )- C /(l + c)]' [6) 
which yields c w 5, and r 2 oo ~ 2 /i _1 Mpc. From this they concluded that there is no conflict 
between the observations of this cluster and the predictions of the standard CDM model. 



3. Consequences for the Cluster Velocity Dispersion 

3.1. NFW profile 

These mass profiles inferred for the cluster CL 0024+1654 based upon lensing 
observations have implications for the velocity dispersion of its dark matter and galaxies if 
the cluster is assumed to be in virial equilibrium. We begin by considering the NFW profile. 
The same N-body simulations of cluster formation in the CDM universe which indicate 
that clusters are, indeed, in an approximate virial equilibrium which can be described by 
the universal mass profile suggested by NFW also yield information about the 1-D velocity 
dispersion ay and its radial dependence. Over a wide range of radii, the halos obtained in 
N-body simulations are roughly isothermal (Tormen, Bouchet, & White 1997; Eke, Navarro, 



-7- 



& Frenk 1998). For comparison with the velocity dispersion of a cluster like CL 0024+1654 
observed within some radius, it is necessary to consider the average velocity dispersion 
of the NFW halo within a sphere of the same radius. The average velocity dispersion of 
CL 0024+1654 was measured by Dressier et al. (1999) to be ay = 1150 km s -1 within 
a radius r ~ 600h~ 1 kpc ps 6r arcs , based upon 107 galaxy redshifts, to an uncertainty of 
roughly less than ±100 km s^ 1 . For the NFW profile proposed by Broadhurst et al. (2000) 
for this cluster, this radius corresponds to r ~ r 2 oo/3. It is a relatively simple matter to 
estimate the predicted velocity dispersion o"v,nfw f° r a given NFW density profile in terms 
of the circular velocity profile f c ,NFW of the same halo, as follows. 

The circular velocity profile of the NFW halo is given by 

2 , , GM(<r) 2 [\n(x + l)-x/(x + l)} , A . 

r x 

where M(< r) is the mass enclosed by radius r, and x = r/rs- The maximum value of 
^c.nfw is f m ax,NFW ~ 0.465 (AttG psrg) 1 / 2 , which occurs at x — x max ps 2.163. The NFW 
velocity profile for the parameters reported by Broadhurst et al. (2000) is shown in 
Figure 1. We obtain f max ,NFW ~ 3340kms _1 for this cluster. According to the detailed 
analysis of numerical N-body results by Tormen, Bouchet, & White (1997), the average 
1-D velocity dispersion within r w r 20 o/3 for simulated clusters, which are well- fit by 
the NFW profile, is somewhat lower than f max ,NFW) but never by more than a factor of 

1.5. This factor of 1.5 agrees very well with the aperture-averaged, line-of-sight ay which 
results from solving the Jeans equation for the variation of the radial ay, r with r inside the 
NFW profile, including the possible effects of anisotropic velocities (Lokas & Mamon 2000; 
Lokas, private communication). Therefore, the NFW profile proposed for CL 0024+1654 by 
Broadhurst et al. (2000) implies an average velocity dispersion for the cluster within the 
radius r ~ 600h~ 1 kpc ps 6r arcs of ov,nfw > ^max,NFw/l-5 ~ 2230 kms -1 , which is too large 
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by a factor of ~ 2 to be compatible with the measured velocity dispersion quoted above 1 . 

3.2. Nonsingular Isothermal Sphere 

The Tyson et al. (1998) mass model in equation (2) does not uniquely predict the 
velocity dispersion oy implied for the cluster even if it is assumed to be in virial equilibrium. 
For one, it is a fit only to the projected density within the arcs at r arcs ~ 100 h _1 kpc and, 
therefore, does not constrain the external mass profile very well. For another, even if we 
restrict our attention to radii within the arcs and assume virial equilibrium, this is still not 
enough to specify oy and its radial dependence there uniquely. As a result, it is not possible 
to compare the predicted ay for the Tyson et al. (1998) mass model with the observed oy 
as directly as we did for the NFW profile above. However, we can at least demonstrate that 
the mass model is plausibly consistent with the observed oy, as follows. 

The projected mass density in equation (2) is well fit by that obtained from the 
truncated, nonsingular, isothermal sphere ("TIS") profile of Shapiro, Iliev, & Raga 
(1999). We discuss the details of the fitting procedure and the results elsewhere (Iliev & 
Shapiro 2000a,b). This TIS represents a particular solution of the Lane-Emden equation 
which corresponds to the outcome of the collapse and virialization of a top-hat density 
perturbation. The size and virial temperature which result are unique functions of the mass 
and redshift of formation of the object for a given background universe. According to this 

1 This assumes p s ~ 0.006 h 2 M Q /pc 3 , the value if p crit (z) is for an Einstein-de Sitter uni- 
verse. Our conclusion that <Jv,nfw substantially exceeds the observed oy remains unchanged 
if, instead, we consider other background cosmologies. The values of tWx.NFW and ov,nfw 
for a low-density universe {VLq = 0.3) are only slightly lower than the values quoted above, 
by 10% for the open case (A = 0) and by 25% for the flat (A ^ 0) case. 
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solution, the central density po is roughly more than 500 times the density at the surface, 
and the core radius r is about 1/30 of the total size. [Note: Our definition of the core radius 
is ?"o,tis = r Kin g /3, where r King is the "King radius" defined in Binney & Tremaine (1987), p. 
228.] As described elsewhere, this solution is a convenient analytical approximation for the 
halos which form from more realistic initial conditions in the CDM model, which reproduces 
many of the average structural properties of the halos found in CDM simulations, except in 
the very inner profile where the TIS model has a uniform-density core instead of a central 
cusp. As such, a fit of the TIS profile to the Tyson et al. (1998) mass model provides a 
plausible, physically-motivated connection between this mass model with a flat density core 
and the implied cluster velocity dispersion. The projected density profile of a TIS sphere 
with a central density p xs 0.064 h 2 M /pc 3 and a core radius r xs 20/i" 1 kpc provides a 
very close match to the Tyson et al. (1998) result for S(r) discussed above. Based on this 
best-fit TIS, the velocity dispersion is 

a V TIS = (AnGporl) 1/2 « 1200 km s" 1 , (5) 

in close agreement with the measured value. 

4. Conclusion 

Two different attempts to invert the observations of strong lensing by cluster 
CL 0024+1654 to solve for the mass profile of the cluster have reached opposite conclusions 
regarding the presence of a uniform- density core versus the acceptability of a central cusp 
like that of the NFW profile. This suggests that, either there is some error in one or both 
of these analyses or else there is some degeneracy in the inversion process which prevents 
a clean discrimination between these different solutions. While we do not claim to address 
the accuracy of either of the two analyses of this cluster, we point out, instead, that there 
may be additional constraints on the allowed mass models which can aid in distinguishing 
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them. In particular, the measured cluster velocity dispersion, which also reflects the mass 
distribution of the cluster, should also be required to be consistent with the mass model 
derived from the inversion of the lensing data. We have shown that this is not the case for 
the cuspy NFW profile which Broadhurst et al. (2000) report is a good fit to the lensing 
data. That profile, if it actually corresponds to a halo formed in the standard CDM model, 
predicts much too high an average velocity dispersion to be consistent with the observed 
value. By contrast, the mass model of Tyson et al. (1998) does not have this problem. 
This suggests that the lensing data for CL 0024+1654 still favor a flattening of the density 
profile at small radii, in conflict with the prediction of cuspy halos by N-body simulations 
of the standard CDM model. 

In view of the importance of this conclusion for the ongoing debate regarding 
the validity of the standard CDM model, it would be valuable if uncertainties in the 
gravitational lense models, such as that due to possible departures from spherical symmetry 
or substructure, were properly quantified. This uncertainty would have to be rather 
extreme, however, to reconcile the Broadhurst et al. (2000) profile fit with the observed 
velocity dispersion of CL 0024+1654. The lensing analysis leads to NFW parameters 
roughly by measuring the mass interior to the arcs, M(< r arcs ) = 47rpsr|/(x arcs ), where 
f(x) = ln(l + x) — x(l + rc)~ 1 and x = r/rs, and determining r$ by matching the logarithmic 
slope 7 of the density profile at r arcs as derived from the projected mass distribution, 
according to :r arcs = — (1 + 7)/(3 + 7). In order to adjust f max ,NFW downward relative 
to the value discussed above by a factor large enough to reconcile the NFW profile with 
the observed velocity dispersion while leaving the measured mass M(< r arcs ) unchanged, 
therefore, r s must be reduced by factor of order 4, to r s 100 /i _1 kpc. This requires an 
observational uncertainty so large as to change the value 7 ~ —1.3 reported by Broadhurst 
et al. (2000) into 7 = —2, which seems very unlikely. 
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Fig. 1. — Circular velocity profiles for cluster CL 0024+1654 implied by the mass profiles 
inferred from the modelling of strong gravitational lensing data by Broadhurst et al. (2000) 
(solid curve labelled "f c ,NFw") an d Tyson et al. (1998) [when the latter is fit by the 
truncated, nonsingular isothermal sphere model of Shapiro et al. (1999)] (solid curve labelled 
"f c ,Tis"-) Horizontal dashed lines show the average velocity dispersions (cry), predicted for 
each mass model, as labelled, along with the value observed by Dressier et al. (1999). Vertical 
dashed line indicates the radius of the arcs, within which the lensing data constrains the mass 
profile. 
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